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Abstract 

Models for topological quantum computation are based on braiding and fusing anyons (quasi- 
particles of fractional statistics) in (2+l)-D. The anyons that can exist in a physical theory are 
determined by the symmetry group of the Hamiltonian. Specifically, any theory of anyons must have 
braiding and fusion rules that satisfy consistency conditions known as the pentagon and hexagon 
equations (collectively, the Seiberg-Moore Polynomial Equations.) Maclane's coherence theorem 
states these are in fact all that is required in order to achieve commutativity of all combinations of 
fusion and braiding (i.e. a consistent physical theory). Two applications of the Hexagon Equation 
yield the Yang-Baxter Equation(YBE) familiar from statistical mechanics: ajaj^iaj = aj^iajaj+i 
where the crj are the abstract braid group generators. It is an unsolved mathematical problem 
to determine in general all the matrix solutions to the YBE. In the case that the Hamiltonian 
undergoes spontaneous symmetry breaking of the full symmetry group G to a finite residual gauge 
group H, however, solutions are given by representations of the quantum double D{H) of the 
subgroup. The quasi-triangular Hopf Algebra D{H) is obtained from Drinfeld's quantum double 
construction applied to the algebra F[H) of functions on the finite group H. As a vector space, 
D{H) = F{H) (g) C[H] = C{H X H) where C[H] is the group algebra over the complex numbers 
and C{H X H) is the space of C- valued functions on H x H. 

A major new contribution of this work is a program written in MAGMA to compute the particles 
(and their properties - including spin) that can exist in a system with an arbitrary finite residual 
gauge group, in addition to the braiding and fusion rules for those particles. We compute explicitly 
the fusion rules for two non-abelian group doubles suggested for universal quantum computation: 
Ss and Ar,, and discover some interesting results, subsystems, and symmetries in the tables. SO{3)4 
(the restriction of Chern-Simons theory S'f7(2)4) and its mirror image arc discovered as 3-particle 
subsystems in the 53 quantum double. The tables demonstrate that both S3 and A5 anyons are all 
Majorana, but this is not the case for all finite groups. (In the appendices, a few other non-abelian 
groups of interest - ^4, A4, and D4 - are included). In addition, the probabilities of obtaining 
specific fusion products in quantum computation applications are determined and programmed in 
MAGMA. Throughout, connections to possible experiments are mentioned. 

Keywords: topological quantum computation, quantum physics, condensed matter 
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I. INTRODUCTION 



Richard Feynman in 1982 was the first to conceive of the idea of a computer based on 
quantum mechanics as potentially offering more computational power than computers based 
on just classical physics. [1] As computer components are beginning to attain distances on 
which quantum mechanical effects can no longer be ignored, continued size reduction would 
necessitate the design of a system that thrives on quantum mechanical distance scales. In 
addition, subsequent developments in algorithms have shown that certain problems believed 
to be classically computationally intractable would be efficiently solvable on a quantum com- 
puter. (Notably, quantum algorithms for factoring large numbers and finding discrete logs 
efficiently |2] would be able to decipher RSA and Diffie-Hellman, a fact which would have 
radical implications for the security of data ranging from credit cards to electronic commu- 
nication). These algorithms would suggest even the asymptotic computational power of a 
device is inextricably linked to the underlying physical mechanisms governing its operation 
in this more general framework. (The Church- Turing thesis asserts that all classical devices 
lead to the same asymptotic complexity.) 

The main idea of topological quantum computation is to exploit in an automatically 
scalable way the physics of the underlying system to obtain computations that are inherently 
protected from errors. Thus the issue of quantum error correction on the software level is 
sidestepped to a certain extent as error-prevention is a physical property of the system 
(effectively hardwired). From the computational perspective, the Threshold Theorem (as 
discovered by A. Steane in [3| and D. Aharonov and M. Ben-Or in |1]) guarantees that 
errorless quantum computation can be achieved using quantum error correction as long as the 
original error rate is below a fixed threshold. However, even if the current threshold estimate 
(1 part in 10~^) is attainable, the overhead required to accomplish the error correction 
procedure is impractical - in general, thousands of qubits are required in order to obtain 
the desired control over one. As it stands, topological accuracy may even be necessary to 
achieve the threshold. 

Information is stored in the global topology of the system and processed via fusion and 
braiding. Since the environment interacts locally and the global topology of the system is 
not effected by local perturbations, the system is inherently protected against decoherence. 
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II. (2+l)-DIMENSIONAL TOPOLOGICAL GAUGE THEORIES SET THE 
STAGE FOR QUANTUM COMPUTATION 



The systems in which topology acts to our advantage are those that are effectively (2+1)- 
dimensional. The striking accuracy that can be achieved in topological systems for even 
classical electronics applications is very promising for the potential accuracy of a topological 
quantum computer. For example, the resistivity in the IQHE is accurate to 1 part in 10^ [5] 
- It was even selected by NIST as the standard of measurement for resistance. Replicating 
that type of inherent accuracy in the actual computations of a quantum computer would 
be an immense advantage over standard classical computations, regardless of asymptotic 
computational complexity. 

In three spatial dimensions, particles can be classified as bosons or fermions according 
to whether they have integer or half-integer spin respectively - those are the only two pos- 
sibilities. However, F. Wilczek showed that in two dimensions the structure can be much 
richer and in fact quasiparticles (anyons) of more general spin and exchange statistics can 
exist. [6] 

The basic mechanism behind Topological Quantum Computation is to interchange and 
fuse these anyons in a controlled manner to achieve computations. The term any-ons, orig- 
inally coined by F. Wilczek, refers to the fact that these quasiparticles are a generalization 
of bosons and fermions in which the particles in abelian theories acquire "any" (rational) 
phase 9 when they are interchanged. Bosons correspond to 9 = Q and fermions correspond 
to 9 = -K. Stability restricts 9 to be any rational value. In non-abelian theories the fields can 
even transform according to a higher- dimensional representation. It is in fact non-abelian 
theories hold the most promise for quantum computation. Anyons have been observed in 
a few different contexts in nature, notably in observations of the Fractional Quantum Hall 
Effect (FQHE) and are postulated to appear in a number of other different models which 
have yet to be experimentally tested. 

Any physical topological theory is specified by the braiding rules and fusion rules. From 
the mathematical perspective, the braiding and fusion rules must satisfy certain consistency 
conditions known as the Seiberg-Moore Polynomial Equations (also known as the pentagon 
and hexagon equations). Maclane's coherence theorem states these are in fact all that is 
required in order to achieve commutativity of all combinations of fusion and braiding (i.e. 
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a consistent physical theory). The two main operations required for topological quantum 
computation are these knotting and fusion rules for the anyons. 

A. Knots 

Knot theory originally began with the suggestion of Lord Kelvin that the fundamental 
building blocks of matter were knots in the ether. While the original incarnation of Kelvin's 
idea was not quite physically grounded, it spawned a research program in mathematics 
that led to many exciting discoveries. As a result, knots have been extensively studied and 
tabulated. The effort received a boost with computer algebra in the 1990s with, among 
others, J. Weeks's computer program to compute various knot invariants. [7] 

With the advent of relativity, physicists showed that space and time are linked and thus 
spacetime in everyday experience is in fact inherently 4- dimensional. However, a modifica- 
tion of Kelvin's idea may have an element of truth in that in (2+l)-dimensions a configura- 
tion of fundamental excitations/particles in particular systems can be described by a knot 
rather than the excitations/particles themselves. Specifically, knotting is one of the basic 
mechanisms by which fundamental excitations can be used to achieve topological quantum 
computation. 

Specializing to particular systems may on the surface appear less fundamental than con- 
sidering the elementary particle constituents of matter in a vacuum. However, it may have 
more to do with particle physics than it would seem on first inspection - In recent years, it has 
been shown that classical notions of the vacuum are false. The vacuum is not a stationary 
medium characterized by the absence of particles and interaction, but rather it is permeated 
by the quark-antiquark QCD condensate and seething with virtual particle-antiparticle pair 
production. In addition, it has been shown that the equations of electromagnetism can be 
obtained as an effective field theory, analogous to those in condensed matter. Thus effective 
field theories may in fact be the only ultimate physical reality. In any event, effective field 
theories play a key role in topological physics. 
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B. Fusion 



Fusion rules are necessary to determine the experimental feasibility of any anyon proposal. 
In particular, if different anyon configurations in the proposal have different cumulative 
long-range fusion rules, the environment will be able to distinguish the states, resulting in 
decoherence. In addition, fusion is the primary means of readout in Kitaev's, Preskill's, 
and Mochon's suggestions for universal quantum computation based on finite groups. Thus 
it is important to understand the fusion rules well in order to build anyon computers. 
Furthermore, the anyon properties of spin, charge, and magnetic flux are determined by 
their particle types (superselection sectors) and these similarly are important to keep in 
mind in the design of experimental realizations of anyons. To this end, this paper will 
culminate with an original program in MAGMA to determine the particles permitted by 
the theory, the properties of those particles, and the fusion rules given the arbitrary (finite) 
symmetry group of the Hamiltonian. 

First, however, it is necessary to understand the theoretical background. Thus the orga- 
nization will be as follows - The physics underlying discete gauge theories will be described, 
including a motivational example of the Aharonov-Bohm Effect in quantum electrodynamics 
(QED) and the mechanism behind how subgroups of the (truncated) braid group B{n,m) 
{n,m G Z>i) of anyons give rise to particles with different charges and fluxes in 2D. The ap- 
propriate mathematical language for discussing braid groups turns out to be quasi-triangular 
Hopf algebras in the quasitensor category structure. The particles are determined by the ir- 
reducible representations of a particular quasi-triangular Hopf algebra obtained by Drinfeld's 
Quantum Double Construction. 

In Section C, a couple of key special cases from the perspective of universal quantum 
computation (6*3 and A^) will be analyzed in depth. The particles (defined by charges and 
magnetic fluxes) and the fusion rules (from the comultiplication of irreducible representations 
of the quantum double) that can exist in these theories will be determined explicitly. In 
addition, the spins of the particles will be computed. A number of interesting results are 
found from the resulting tables. In particular, some subsystems are discovered that are 
interesting from either a mathematical or quantum computation perspective. The emphasis 
is on theories that are physically relevant. The actual MAGMA code and the results for 
more examples are provided in the Appendices. An overview of the MAGMA program and 
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a summary of the main results is provided. 

Finally, possible experimental implementations and applications to universal quantum 
computation will be discussed. The probabilities of obtaining certain fusion products in the 
quantum computation scheme will be derived. The resulting probabilities are implemented 
in the MAGMA code included in Appendix A. 

III. MOTIVATION: TOPOLOGICAL EFFECTS IN PHYSICS 

Conservation laws are critical to many physics analyses as they often simplify calculations 
dramatically and aid intuition. Noether's Theorem states that for every symmetry of the 
Hamiltonian, there is a conserved current and provides an explicit formula for calculating the 
conserved current. However, it turns out that this list is not exhaustive - In fact, topological 
conservation laws can arise that cannot be traced to an underlying Noether current. 

There are two major manifestations of topological effects - in coordinate space, the 
Aharonov-Bohm Effect, and in momentum space, the Berry phase. 

A. Aharonov-Bohm Effect 

As motivation, recall the Aharonov-Bohm effect in quantum electrodynamics. [8] Initially, 
gauge fields were introduced into classical physics as auxiliary fields to aid calculations - 
While they were present in the canonical formalism, they were not necessary and did not 
appear in the classical equations of motion (Maxwell's equations in the case of electromag- 
netism). However, quantum mechanics promotes observables to operators and mandates 
the use of the canonical formalism. As a result, the topology of the gauge fields can indeed 
have observable physical consequences. In particular, in quantum mechanics, electrons be- 
come waves which exhibit an interference pattern. Aharonov and Bohm predicted that a 
change in the gauge fields (by for example, the introduction of a solenoid in between the two 
beams) will cause a shift in the interference pattern observed. As quantum mechanics only 
permits local interactions among fields and the two electron beams do not pass through the 
region in which the gauge fields are nontrivial, the use of the vector potential in this case is 
inescapable. 

While first introduced as a Gedanken-experiment for demonstrating the infiuence of topo- 
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logical defects on physical systems, the Aharonov-Bohm Effect has now been verified by 
many actual physical experiments (including a 1998 experiment using magnetostatic metal 
rings [9]). The typical experiment demonstrating the Aharonov-Bohm effect is the 2 slit 
electron diffraction experiment in which electrons are fired at the slits one at a time and 
the diffraction pattern on a screen a fixed distance away is measured. A change in the 
interference pattern is observed when a solenoid (ideally, infinitely long so the magnetic 
field completely vanishes outside the solenoid) is placed in between the two slits. Thus it is 
found that the vector potential itself has intrinsic physical meaning in this effectively two 
dimensional setting, not just the curl as originally supposed in classical electrodynamics. 

To formalize the exact effect, apply Feynman's path integral method to the system at 
hand. The phase for the interference term is determined by the path difference between the 
two electron beams. The degree of the effect is thus determined by the integral of the vector 
potential around a loop surrounding the solenoid. 

Given a non-simply connected system, if the gauge potential is changed by an amount 
6A, then the resulting change in the wave equation is only invariant if f 6A{r) ■ dr = n-^ 
where n is an integer. Thus the fractional part of § A{r) ■ dr has measurable physical 
consequences. 

In broken gauge theories, the flux is quantized - In order to ensure the wave function is 
single valued, the vector potential can only encircle the solenoid an integer number of times. 
Note that this is a purely topological effect as it is determined solely by the winding number 
of the path. (In mathematical terms: the path integral exponent corresponds to an element 
in the fundamental group of the space . Specifically, the solenoid is a puncture in the plane 
(the third dimension is trivialized) so 7ri(]R^ \ x M) = 7ri(S'^) = Z.) Thus topology can play 
an influential role in determining quantum physics phenomena. In fact, since topological 
interactions are by nature robust against local perturbations, harnessing this resource for 
quantum computation could lead to a fault-tolerant quantum computer. 

B. Flux Metamorphosis 

In this paper, the topological effects arising from Hamiltonians in (2+l)-dimensions with 
finite symmetry groups will be considered. Let H be a nonabelian group. The basic idea 
is as follows: Consider two particles with fluxes in the residual gauge group: hi,h2 € H. 
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FIG. 1: Two particles undergo flux metamorphosis when exchanged. All nontrivial gauge trans- 
formations occur in the "Dirac strings" of the particles, denoted here by the lighter lines going off 
to infinity, (a.) /ii,/i2 £ H (b.)-(c.) /12 is moved counterclockwise around hi. (d) = hih2h^^ ■ 

Fluxes introduce a nontrivial holonomy on the set of ground states that effects any test 
particle whose path winds around the flux. WLOG take all the nontrivial parallel transport 
to occur in a single strand, the Dirac string attached to the particle. (This corresponds to 
fixing the gauge.) The Dirac string will be represented as a single thick line in the vertical 
direction for all particles. 

Now, the long range flux of the two particles combined is hih2. Exchanging the two 
particles counterclockwise, the particle with flux /12 passes through the Dirac string of the 
particle with flux hi. Thus the flux of the particle with flux h2 will change to a new flux, call 
it /i2- Interchanging the particles should not have an effect on their cumulative long-range 
topological quantum numbers. Thus, when the two particles are exchanged counterclockwise, 
the cumulative flux of the two particles should remain equal to hih2. Thus the relation 
/ig/ii = hih2 is obtained, or /ig = hih2h^^. This effect is called flux metamorphosis and 
applies to any 2D topological theory based on a non-Abelian group. 

IV. PHYSICS: HAMILTONIAN DESCRIPTION 

Modern physics has evolved closer to group theory as Hamiltonians can be classified solely 
on the basis of their symmetry groups and particle content. In fact, the current standard 
model is based on the Hamiltonian described by the gauge symmetry group SU (3) x SU (2) x 
f/(l). The symmetry groups seen in nature are given by spontaneous symmetry breaking 
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(SSB) of this gauge group to subgroups. SU{3) is the strong force and SU{2) x U{1) is 
the electroweak force. Electromagnetism arises in SSB of the SU{2) x C/(l) part to a C/(l) 
subgroup which involves both SU{2) and f/(l) components. (The bosons of the weak 
force are charged, so the U{1) cannot be cleanly separated from the SU{2) to be interpreted 
as charge.) In fact, superconductivity (also likely to be critical for quantum computation) 
corresponds to the further SSB of the U{1) electromagnetic gauge group to a subgroup. 
This section will be devoted to understanding the underlying Physical foundations behind 
the existence of the particles. First, SSB will be defined and a generic Hamiltonian for 
a physical system is provided. Then the restriction is made to the case of n particles in 
(2 + l)-dimensional spacctimc and the symmetry groups that can occur for this restricted 
class of systems are determined. 

A. Hamiltonian 

Since the Hamiltonian of the system as a whole has to transform as a singlet under the 
representations of the symmetry group, the representations limit the types of terms that can 
occur in the Hamiltonian. In fact, the specification of a symmetry group and the particles 
that actually are present determines the Hamiltonian uniquely. In (2+l)-dimensional gauge 
theories based on finite groups, the particles must satisfy consistency conditions so that 
braiding and fusion between three particles commute. The consistency requirements restrict 
the types of particles that can possibly exist in the theory. 

To be explicit, the generic action in the Yang- Mills- Higgs theory can be written as 

S — SyMH + Sjnatter 

In quantum field theory there are two types of symmetries: 

1. global symmetries (physical symmetries of the theory, which have an associated con- 
served Noether current) and 

2. local gauge symmetries (nonphysical redundancies in the coordinates chosen, so do 
not have a conserved Noether current). 

Suppose the action S (equivalently Hamiltonian, Lagrangian) has symmetry group G. 
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In the former case 1, if the global symmetry is broken to a finite subgroup then according 
to Goldstone's theorem a massless particle (a Goldstone boson) is created for every broken 
symmetry generator. The emphasis in this thesis is on the latter case 2, in which the vector 
gauge fields (usually originally massless) acquire mass through the Higgs mechanism. 

B. Symmetry Breaking 

Starting with a Hamiltonian with gauge group G, there are two basic ways to break 
symmetry: 

1. Add an external potential: If the potential does not respect the full symmetry of the 
original Lagrangian, then the symmetry of the resulting Lagrangian becomes some 
subgroup of the original symmetry group. 

2. Spontaneous Symmetry Breaking (SSB): (described in most Quantum Field Theory 
textbooks, see for instance [TU]) - In this case, the full symmetry of the original La- 
grangian is broken by the choice of a particular ground state (frequently this is referred 
to as the vacuum expectation value, v.e.v.) The order parameter space is isomorphic 
to G/H (an element of G/H rotates from one ground state to another in the system.) 

In the Standard Model of particle physics, the gauge boson acquires mass through the 
Higgs mechanism. The same phenomenon occurs in Condensed Matter Physics in supercon- 
ductivity - in fact, an (equivalent) formulation of superconductivity is that the photon field 
of electromagnetism acquires mass within the superconductor. 

C. Monopoles and Instantons 

Monopoles can arise as either Dirac or Polyakov-'t Hooft monopoles in field theories. 
Both carry quantized flux, yet there are some crucial differences to keep in mind. ( is a 
good reference.) 

1. Dirac Monopoles 

P. Dirac predicted the existence of monopoles in the quantum theory of electromagnetism 
as a method of accounting for charge quantization. These Dirac monopoles are singular 
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solutions of the field equations carrying infinite energy. In addition, Dirac monopoles have 
the aesthetic appeal of making Maxwell's equations symmetric. To date, however, magnetic 
monopoles have not been observed in nature (Every bar magnet appears to split into smaller 
bar magnets when torn apart down to the smallest scales measured). Dirac monopoles 
require the introduction of non-electromagnetic degrees of freedom, charged particles whose 
mass and spin are free parameters. 

Mathematically, the monopoles are labeled by elements of the fundamental group 
MG/H)). 

2. Polyakov- 't Hooft Monopoles 

In broken gauge theories, it is not Dirac monopoles but rather Polyakov-'t Hooft 
monopoles that play a key role. These are of a very different nature than Dirac monopoles 
- they are non-singular solutions of the field equations that carry finite energy, and neces- 
sitate no new degrees of freedom, i.e. all their parameters can be expressed in terms of the 
ordinary scalar and vector meson parameters. 

Recall that TTi{M),i > 2 is always abelian. 

3. Instantons 

When the gauge group G is not simply connected, the lift H of the subgroup H to the 
universal cover determines the different magnetic vortices. 

D. Higgs Screening 

In light of the preceding discussion of the Aharonov-Bohm effect, it is important to note 
that the physical charge de facto plays two distinct roles in the Yang Mills Higgs theory. - 
It serves as both the coupling constant for the Coulomb interactions and as the coupling 
constant for the Aharonov-Bohm scattering. In the Higgs phase. Coulomb charges are 
screened at large distances (compared to the scale where Ma is the mass of the Higgs 
gauge fields). However, topological interactions may persist at those distance scales. [12] 
In order for a given charge q to be detected via Aharonov-Bohm scattering, the screening 
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charge g' must not couple to the Aharonov-Bohm interaction. 

The resulting massive vector bosons mediate short-range interactions only (the field 
strength decays exponentially as the mass of the Higgs field). Since we are interested in 
the global properties of the system, consider the long-distance (low energy) regime in which 
the Higgs field is condensed in one of its ground states. Since the field strength decays ex- 
ponentially, the interactions between two particles due to their fundamental gauge quantum 
numbers (throughout this paper the term electric charge will be used, by analogy with the 
electromagnetic case) at large separation are exponentially suppressed. (In the case of the 
electromagnetic field, this can be interpreted as Coulomb screening by the Higgs medium.) 

However, long range topological quantum numbers (which will henceforth be called flux) 
are still well-deflned and the system can thus be manipulated by the braiding and fusion 
determined by these topological quantum numbers. Since the resulting operations depend 
only on the topological parameters of the particles, the specific process or path by which 
the braiding and fusion is carried out does not effect the conclusions. (The framework 
described can thus be applied to many different possible topological quantum computation 
implementations. Since the computation is robust against local deviations, this would also 
lead to reliable quantum gates.) 

E. Classification Based on the Symmetry Group 

1. Three Dimensions 

Consider a physical system consisting of n identical point particles. (Note that for these 
particles to be truly indistinguishable, the system has to be quantum mechanical.) In three 
or more dimensions, any loop that circles twice around a point can be continuously deformed 
to the identity. (7ri(S'0(3)) = Z/2Z) Thus the only physical property corresponds to the 
particular permutation of the n particles (the winding number is always trivial), and the 
fundamental group of the configuration space of the underlying manifold of this physical 
system is given by 7ri(C„(M)) = 7ri((M" — D)/Sn) ~ S'„, where D denotes the configura- 
tions in which two or more particles coincide. Thus the Hamiltonian transforms under a 
representation of the symmetric group. The symmetric group Sn always has two unitary 
irreducible representations (UlR's) of dimension 1 - the trivial representation and the anti- 
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symmetric (sign) representation. In ordinary three dimensional space, these correspond to 
Bose-Einstein and Fermi-Dirac statistics respectively. Fractional statistics could in theory 
arise in three dimensions, but would have to correspond to a higher dimensional represen- 
tation of the symmetric group. However, as explained by Haag in any physical theory 
including fermions obeying parastatistics in three dimensions can be modeled as a theory 
including only regular fermions by including some hidden variable, without changing any of 
the resulting physics. 

2. Two Dimensions 

Indistinguishable Particles: Braid Group On the other hand, in two dimensions 
the fundamental group of the configuration space of a system of n identical particles is given 
by the more general braid group: 

7ri(C„(]R^)) ~ 5„(]R^) = {(Ti, 1 <i <n - l\(Tiai+iai = ai+iaiai+i, 

aiaj = (Tj(Ti\f\i - j\ > 2). 

The (jj are defined physically by interchanging the i^^ and {i + 1)'^* particles counterclock- 
wise. In contrast to the symmetric group, not all braid group generators satisfy af = Id. 
Intuitively, two dimensional space allows for more possibilities because particles are topo- 
logically pointlike. Adding a particle to a 2D surface is like poking a hole in the surface, so 
loops around the hole can no longer be continuously deformed to a single point. (Note that 
our motivational Aharonov-Bohm effect example simulates 2D as the solenoid is infinite in 
length, so only 2 of the dimensions really matter. This corresponds to the fact that the 
loops around a line are characterized by the fundamental group (first homotopy group), 
7ri(M2 \ line) = ni{S^ xR) = Z.) 

Truncated Braid Group The braid group as defined is infinite and torsion-free (all 
elements of the group have infinite order). For theories with a finite symmetry group, it is 
actually the truncated braid group B{m, n), that governs the interaction of particles. These 
are the braid groups on n elements with the additional condition that cr™ = 1 for some 
finite m (m is taken to be minimal, m = 2 corresponds to the case of Sn, so the truncated 
braid group is indeed a generalization of the symmetric groups.) When a acts on physical 
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particles, the Hamiltonian transforms as a representation of a (the braid operator), which 
we call the R matrix. (More on this topic will be discussed in the next section.) 

Distinguishable Particles: Colored (Pure) Braid Group If the particles are 

distinguishable, the configuration space is given by C{{Ry'-D). The colored (pure) braid 
group P(m, n), for n even, describes the fundamental group of the configuration space of the 
underlying manifold. (Intuitively, only monodromy operators make sense since after each 
exchange the two particles have to be back in their original positions.) Formally, the pure 
braid group P{m, n) is the kernel of the homomorphism (p : Bn ^ Sn mapping each braid 
to the associated corresponding permutation. 

Partially Colored Braid Group If the system consists of a combination of distinguish- 
able particles and indistinguishable particles, then the fundamental group is the partially 
colored braid group. 

3. Monodromy 

The monodromy operator is given by two applications of the braid operator, i.e. it 
interchanges the particles twice. Since the particles are now back in their original positions, 
this is equivalent to winding one particle counterclockwise around the other. In the case of 
nonidentical particles, the a operator, and thus the R matrix, is not well-defined as after 
the interchange, the particle in each slot is fundamentally different from that before the 
exchange. However, the monodromy operator still makes sense. Now, if the particles are 
both either bosons or fermions then the action of the monodromy operator is trivial, so in 
some sense this is an indicator of whether or not anyonic statistics are present in the theory. 
In quantum interference experiments and the Aharonov-Bohm scattering experiment, the 
monodromy operator, not the braid operator, is actually measured in the laboratory. It is 
sometimes easier to calculate the monodromy operator as it can be written solely in terms 
of the spins of the particles. 
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V. SPIN-STATISTICS CONNECTION 



A. The Abelian Case 

While the traditional proof of the spin-statistics connection due to Wightman provides a 
mathematically rigorous proof based on the axioms of local relativistic quantum field theory, 
the approach in this paper provides an intuitive topological proof that makes evident the 
fact that the spin-statistics connection is in fact a topological property that depends only 
on the existence of antiparticles and Lorentz invariance. (This proof was originally proposed 
by Finkelstein and Rubinstein in the case of three spatial dimensions, then adapted to 
the two dimensional case by Propitius de Wild. [14J) Consider the process in which two 
particle/antiparticle pairs are created from the vacuum, the particles are interchanged, then 
the new pairs annihilate. The associated spacetime ribbon can be deformed via Lorentz 
transformations into the ribbon for rotation of a single one of the particles by an angle of 
2tt about its center, as can be seen from diagram [2]. Thus from the topological perspective, 
these two processes are seen to be physically equivalent. 

Thus the statistical phase factor e*^ picked up by exchanging the two particles must in 
fact be equal to the spin factor e^'^** obtained by rotating a particle of spin s. Note that the 
procedure in the above proof is independent of whether the particles or the antiparticles are 
interchanged (the same orientation results either way). Thus as an additional consequence 
we obtain that the spin of the antiparticle is equal to that of the corresponding particle. 

B. The Non-Abelian Case: Generahzed Spin-Statistics Connection 

When the group is non-abelian, the canonical spin statistics connection must be general- 
ized as the topological equivalence in Fig. [? ] shows. Specifically, when a composite particle 
C (the result of fusion of A and B) is rotated, not only is the phase of the composite picked 
up, but each of the constituents A and B are also rotated. Thus both A and B must both 
be rotated back in the opposite direction, each picking up a phase. 

jy-ABC jd2 _ „2TTi{s^c,j)-S(A,a)-SB.p) L^ABC 

When the two particles are the same, a square root version of the generalized Spin- 
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FIG. 2: Note that the effect of interchanging two particles in two separate particle/antiparticle 
pairs can be continuously deformed into the process under which one of the particles is rotated by 
an angle of 27r about its center. Thus the two effects should be topologically the same. 

Statistics Connection is satisfied: 

1{-AAC T3 _ ^p7rj(s(c7,T,)-2s(^_Q)) 7>-AylC 

VI. DETERMINING THE PARTICLES FROM THE HAMILTONIAN 
A. Abelian 

Consider the following action for a minimally coupled scalar field with an Abelian sym- 
metry group: 

S = SyMH + Smatter = J x{-^F^,F'^'' + (D^$)*D^<I> - - jM^) 

where the covariant derivative is given by -Dp$ = [dp + iNeAp)^ (i.e. the Higgs field has 
charge A^e with respect to the compact U{1) symmetry). The last term represents the 
contribution from the matter fields while the others are the Yang- Mills- Higgs contribution. 
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FIG. 3: The Generalized Spin-Statistics Connection applies to two particle fusion in non-abelian 
theories. On the left, fusion is performed after a monodromy operation. This can be continuously 
deformed via Lorentz transformations to the diagram on the right, in which the particles are each 
twisted by a — 27r rotation before being fused, then the fusion product is twisted by 27r. 

Terms in the Lagrangian density that can be written as a total derivative are not included 
as usually the terms cancel at the boundaries. (This is not always the case. In particular, 
Chern-Simons theories contain an additional term which is a total derivative but the bound- 
ary conditions still contribute precisely because there is a non-trivial topological winding 
at the boundary as in the case, for instance, of axions. There is also a description of the 
Fractional Quantum Hall effect in terms of Chern-Simons theories. However, in this paper, 
the emphasis is on topological effects due solely to the ordinary Lagrangian terms.) 




B. Non-Abelian Cases 



In the case of the nonabelian Higgs field, the field tensor must be given an extra index a. 
Symh = j ^^^(-^F'^^^F;, + (D^<l>)t . D,^ - 
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where the covariant derivative is given, for instance, in the adjoint representation of 
the nonabehan group by = 9^0^ + gf°^^'^A''^^c- Exphcitly, the field strength is 

C. Determining the Particles 

Note that there is spin-charge separation in the system as the flux and the charge can 
be measured independently (and there are both charged particles with no magnetic flux 
and magnetic fluxes with no electric charge). A specific particle type is determined by a 
conjugacy class A of the symmetry group H and a representation a of the centralizer of a 
representative of that conjugacy class. (The centralizers of each element of the conjugacy 
class arc isomorphic, so the representation thus defined is independent of the particular 
representative chosen.) In particular, label the elements of each conjugacy class "^C — 
{^hi,^ /i2, ■■■i^ hk}. Let "^N denote the subgroup of H that is the centralizer for ^/ii. Choose 
a set of representatives of H/^N — {""vi, Vk}- The states 

{\''h,rv,)} 

then constitute a basis for the internal Hilbert space V^. 

A bit more formally, given the fiux values take values in a group H, then quantum 
mechanically we should consider complex combinations C[i7]. The electric charge takes 
values in the dual F{H). The particles thus correspond to representations of Drinfeld's 
quantum double construction applied to the Abelian algebra F{H) of functions on a finite 
group H. Generators are given by combinations {P^,g}h^g^H-i where Ph is the operator 
projecting the fiux h & H and g & H corresponds to the global symmetry transformation. 
The operators Ph are projections as PhPw = Sh,h'Ph- We demonstrated earlier that the 
action of the global symmetry transformations does not commute with the fiux projection 
operators, and in fact: gPh — Pghg-ig due to fiux metamorphosis. Combining these two 
relations we obtain the multiplication rule for two elements of the quantum double: Phg ■ 
Pwg' — ^h,gh'g-^Phgg- The particlcs permitted by a discrete gauge theory are given by unitary 
irreducible representations of the quantum double. Note that since the quantum double is an 
algebra, not a group, representations 11 here are defined not only to respect the multiplicative 
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structure, but are also required to be linear, i.e. addition and scalar multiplication are 
preserved. Each unitary irreducible representation of D{H) is labeled uniquely by fixing a 
conjugacy class of the underlying group H and specifying the representation of the centralizer 
of a representative element in that conjugacy class. (The centralizers of all elements in a 
conjugacy class are isomorphic, so it does not matter which representative element is chosen.) 
The algebra is semisimple - it can be decomposed into a direct sum. Thus the sum of squares 
of internal Hilbert spaces = the dimension of the quantum double. (This will be used as a 
consistency check to confirm the explicit dimensions obtained in future sections.) 



D. Changing Basis: The S-Matrix 

The S-matrix implements the change of basis between the given basis of the fusion algebra 
and the diagonal basis. (Such a matrix S must exist since the algebra is commutative and 
associative.) Explicitly, the S matrix can be written as: 



S^^ -^tr R 
l-ti 



al3 



' ' ^hie^CBhjeBC:[^hiBhj]=e 

In the appendices, a MAGMA program that computes the S matrix values is included. 



E. Spin and the T-Matrix 

The diagonal T- matrix records the information about the spins. Explicitly, the T matrix 
is given by: 



Within each internal Hilbert space, all particles carry the same spin. Furthermore, the 
only particles which have nontrivial spin are dyons with both nontrivial flux and nontrivial 
charge. Note that in the construction for the quantum double, by definition the conjugacy 
class A representative element ^/ii commutes with all elements in the centralizer ^A^, so 
Schur's lemma yields: 
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a 



for every irreducible representation a. This formula is implemented in the MAGMA 
program to yield the values of the spin for the various anyonic particles. While this does 
not distinguish between spin s = and s = 1, it does the best that is possible in any 
superselection sector as it is always possible to add orbital angular momentum. In the 
specific groups discussed in the next section, a number of exotic spins are discovered. 

F. Magnetic-Electric Charge Duality 

The symmetry transformations of the second factor (C(if)) in the quantum double con- 
struction relate D{H) to the electric charge. Similarly, the first part F{H) of D{H) yields 
the magnetic fluxes. Note that these are in fact dual to each other, as F{H) CH*. The 
magnetic-electric charge duality is implemented by the S'L(2,Z) action generated by the S 
and T matrices. The action of the 5* matrix has an interpretation as a Fourier transform. [IS] 



Similar to regular quantum field theory of the collision of two particles in three dimensions, 
the fusion of two particles can lead to multiple possible different outcomes. The fusion rules 
are the analog in this two dimensional setting of the selection rule for addition of angular 
momenta Ji, J2 in quantum mechanics which permits values between —\J2 — Ji\ and J2 + Ji 
as feasible outcomes for the combined system. The fusion coefficients indicate how many 
different ways there are to obtain that outcome. (However, this still does not tell us the 
probabilities of obtaining the specific outcomes, so it is not quite the analog of the Clebsch- 
Gordon coefficients.) 

The fusion coefficients are encoded in the comultiplication of the dual D{H)*. Explicitly, 



1. Fusion 
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where as before the modular S matrix is given by: 



lAB 

'a/3 



1 



.-2AB 
a/3 



H 



tr R 



1 



E 



tr{a{^xr^''hfxi))*tr{P{''xj'''hfx,))* 



'^hi&'^CBhj&BC:[AhiBhj]=e 



2. Antiparticles 

An antiparticle is defined to be the representation which when fused with the original 
particle yields the identity as one possible outcome. It turns out that the antiparticle is 
unique, and the charge conjugation operator C — takes a particle to its antiparticle. 
Note that since the S matrix commutes with the T matrix (which encodes the spin), the 
spin of a particle is the same as that of the antiparticle. 

G. Alice Physics and Cheshire Charges in Non-Abelian Gauge Theories 

Cheshire charge is a non-localizablc form of charge that is carried by a configuration 
of particles. Globally, charge must be conserved. Thus, a special type of entanglement 
is observed if the particles with Cheshire charge are kept far apart from one another so 
that charge cannot be transferred by local interactions. This entanglement opens up the 
possibility of topological quantum computation using anyons. In particular, consider a 
collection of particle/antiparticle pairs. Each local particle/antiparticle pair has trivial total 
fiux and charge - If kept sufficiently separated, the environment cannot distinguish between 
states with long range quantum numbers the same as those given by the ground state. 
However the configuration still has nontrivial Cheshire charge, so the state cannot annihilate. 
In other words, since the environment acts locally on the system, it cannot effect the non- 
local Cheshire charge carried by the configuration. Thus the system is protected against 
decoherence. 



22 



VII. SUMMARY 



The topology of the configuration space restricts the possible symmetry groups of the 
Hamiltonian. In turn, every term in the Hamiltonian must transform as a singlet under 
all the representations of the symmetry group so this restricts the types of terms that 
can appear in the Hamiltonian. In (2+l)-dimensional systems, the truncated braid group 
must be used instead of the symmetric group. Consistent braiding and fusion then impose 
conditions which must be solved by any particles in the theory. Particles have two degrees 
of freedom given by the magnetic flux and electric charge. In (2+l)-dimensional theories 
based on a finite symmetry group, representations of the quantum double of the symmetry 
group determine the kinds of particles that can possibly exist in the given theory. The fusion 
rules are then encoded by the comultiplication of the dual. 

VIII. EXPLICIT DERIVATIONS OF PARTICLES AND FUSION RULES FOR 
SPECIAL PHYSICAL THEORIES 

A. Abelian vs. Non-Abelian: Physical Origin and Computational Potential 

In Abelian groups, the fusion rules obtained by the quantum double formalism are deter- 
ministic as only one product can result from fusion of two particles. However, in non- Abelian 
Groups, the fusion rules obtained for at least two input particles include multiple possible 
output particles so the product is not determined solely by the symmetry group. (Letting 
-^AB the number of ways of obtaining particle C when the particles A and B are fused, 
mathematically this condition is just YliC ^ab > 1 some particles A, B.) To obtain a 
deterministic result, the more detailed dynamics of the specific interaction have to be taken 
into account. 

Prom another perspective, non- Abelian groups lead to multidimensional Hilbert Spaces so 
there is more wiggle room for changes in the state (i.e. the operators can be noncommutative 
in a multidimensional space) . This is the key behind the computational power of non- Abelian 
groups in quantum computation. 

Abelian groups are severely limited in their applications to quantum computation by 
the fact that they can only be used to store information (not to process it). However, 
Abelian groups nonetheless underlie some physical systems and are used in Kitaev's Toric 
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Code for storage. [16] Specifically, superconductivity in a material corresponds to SSB of the 
electromagnetic U{1) gauge field - the standard condensed matter Cooper Pair, for example, 
arises from the SSB of the electromagnetic U{1) to Z/A^Z with N = 2. 

B. Applying the Quantum Double Formalism to Abelian Finite Groups of the 
Form Z/NZ 

Since finite Abelian groups are the simplest example of the theory, we will begin there. 
In this section we consider the effect of spontaneous symmetry breaking of an Abelian U{1) 
gauge group down to the finite subgroup Z/NZ. Note that screening precludes electromag- 
netic effects between two particles from influencing the long-distance interactions. However, 
the topological charges persist. For these Abelian theories, the monodromy operator is just 
a phase. 

1. Physical Theory 

The spectrum in Z/NZ theories consists of pure charges, pure fluxes, and dyons produced 
by fusing these charges and fluxes. Let 



denote a flux a G Z/NZ charge n G Z/NZ particle. 
2. Braiding and Fusion 

The braiding and fusion rules for Z/A^Z theories are particularly simple. Braiding two 
identical particles, one encounters the Dirac string of the other, which introduces a phase: 



a)\n) 



a, n) 



R\a, n)\a, n) = e 



2iri 
N 



(na) 



a, n)\a, n) 



For non-identical particles, only the monodromy operator makes sense: 



i?>,n)|a',n') = e'^('^»'+"'") 



a, n) |a', n) 
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The fusion rules for this theory amount to summing the charges and fluxes mod N: 

\a,n) X \a',n') = \ [a + a'], [n + n']) 

Instantons in effect perform the modulo N summation. The charge conjugation operator 
takes each charge, flux combination to the additive inverse charge, flux mod N. This in fact 
corresponds to the antiparticle. Note that the only Z/NZ gauge theories which contain a 
Majorana particle (i.e. a particle that is its own antiparticle) are those for which N is even. 
In that case, the only Majorana particle has both flux and charge y. 

The scattering cross section for two distinguishable particles is given by: 

da sin^ ( {na' + n'a) ) 
d9 " 27rpsin2(|) 

where p is the relative momentum of the two particles, and 9 is the scattering angle. 
Comparing this with our expression for the monodromy operator R^, we have thus verified 
that the cross section only depends on the monodromy operator for the Z/NZ theory. 

In the case of two identical particles, the cross section must be modified to take into 
account quantum exchange: 

da _ sm\^) sm\^) 
d9 ~ 27rpsin2(^/2) ^ 27rpcos2(^/2) 

C. Applying the Quantum Double Formalism to Specific Non-Abelian Finite 
Groups 

1. S3 

The conjugacy classes of the symmetric groups arc determined by the cycle decompo- 
sition. In particular, for 5*3, there are 3 conjugacy classes with representative elements (), 
(12), (123) - These correspond to fluxes in the spontaneously broken gauge theory. The 
centralizers of a representative element of each conjugacy class are given by: S3, Z/2Z, 
and Z/3Z, respectively. The possible charges are given by representations of the centrahzer 
corresponding to the particle's flux. 

The pure charge sector corresponds to flux () and representations of the centralizer ^3. 
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Since S3 has three conjugacy classes, there are three irreducible representations of 5*3, with 
dimensions 1,1, and 2 (the unique unordered solution {rii, 712, 713} to the equation + n| + 
77,3 = 6 = IS'al). These are given by the identity representation, the sign representation, and 
a 2 dimensional representation. 

For 5*3, the particles with nontrivial flux can have charges that are representations of Z/2Z 
and Z/3Z respectively. These last two centralizer groups are Abelian - the corresponding 
representations are thus just roots of unity. 

• Flux: Identity Element (): In this case, the charges are given by representations of 
the centralizer S3: 



Representations of S3 



Repn 





(12) 


(123) 


Xi 


1 


1 


1 


X2 


1 


-1 


1 


X3 


2 





-1 



• Flux: Two-cycle Conjugacy Class (Representative Element: (12)): In this case, the 
charges are given by representations of the Centralizer Z/2Z. 



Representations of Z/2Z 



Repn 





(12) 


Pi 


1 


1 


P2 


1 


-1 



• Flux: Three-cycle Conjugacy Class (Representative Element: (123)): In this case, the 
charges are given by representations of the Centralizer Z/3Z. 



Representations of Z/3Z 



Repn 





(123) 


(132) 


Pi 


1 


1 


1 


P2 


1 


g27ri/3 


g47ri/3 


P3 


1 


g47r«/3 


g27ri/3 



In summary then, there are 8 different types of particles, given by: 
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A 


B 


c 


D 


E 


F 


G 




A 










E 




G 


H. 


& 


B 




' 1 




D 




G 


H 


C 








[HE 


[HE 


G+H 


F+H 


F+G 


D 


D 


E 


D+E 


A+C+F+G+H 


E+C+F+G+H 


D+E 


D+E 


[HE 


E 


E 


D 


D+E 


B+C+F-H3+H 


A+C+F+G+H 


[HE 


[HE 


D+E 


F ^^^^1 


G+H 


CHE 


[HE 




C+H 


C+G 


G 


G 


G 


F+H 


D+E 


&tE 


C+H 


A+B+G 


C+F 






H 


F+G 


D+E 


&tE 


C+G 


C+F 


A+B+H 



FIG. 4: Fusion Rules for ^3 



Particle 


Flux 


Charge 


Spin (e2^^"(-4,-)) 


Repn Dim 


Dim Internal Hilbert Space 


A 





Xi 





1 


1 


B 





X2 





1 


1 


C 





X3 





2 


2 


D 


(12) 


Pi 





1 


3 


E 


(12) 


P2 


1 

2 


1 


3 


F 


(123) 


Pi 





1 


2 


G 


(123) 


P2 


1 

3 


1 


2 


H 


(123) 


P3 


2 
3 


1 


2 



where the Spin s is chosen in the range < s < 1. 

Note that particles E, G, and H all have fractional spin - particle E is a regular fermion, 
but G and H have spin | and | respectively, so this is an explicit demonstration of fractional 
statistics. As a check to this table, we do indeed have the equivalence that the sum of the 
dimensions of the internal Hilbert spaces is 1^ + 1^ + 2^ + 2 * 3^ + 3 * 2^ = 36 = jS'sp, the 
dimension of the quantum double. 

The monodromy operator is given by the square of the braiding operator. 

Fusion Rules 

Note that in the table for 5*3, every particle is its own antiparticle. Also, there are many 
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subsystems which are closed under fusion. In particular, the subsystems 

a,B}, {A,B,C}, a,B,G}, {A,B,F}, a,B,H}, a,B,F,G}, a,B,F,H}, a,B,G,H}, 

{A,B,F,G,H} and {A,B,C,F,G,H} 

are all closed under fusion. Thus if not all the particles were initially present, then they 
need not arise by either fusion or braiding. However, as soon as D or E particles are present 
in the system, with some probability all of the other particles will appear. In addition, 
quantum mechanically virtual pair production can occur so all the particles appear with 
some probability, if only for a short time, regardless of the initial particle content. 

Note that the {A,B,C} subsystem may be physically natural in a theory in which the 
pure electric charges A, B and C are much lighter than the rest of the particles so we can 
integrate the others out. This is potentially reasonable because magnetic fluxes tend to 
be much heavier than their trivial flux counterparts. The three-particle subsystems will be 
analyzed further in the section on quantum computation. 

Finally, there are a number of unexplained symmetries in the fusion table that a priori 
would not necessarily be expected based on the 5*3 symmetry group of the underlying theory. 
In particular, the fusion table is invariant under any permutation of the labels C, F, G, and 
H. Physically, however, the C particle tends to have a much smaller mass than the F, G, 
and H particles since it is a pure electric charge. 

5*3 Symmetry Group Implementations 

Note that S3 is isomorphic to Cs^, the symmetry group of a honeycomb lattice as observed 
naturally in, for instance, graphene. Thus it is perhaps reasonable to expect to find an 
^3 system conducive to topological gauge theories in nature. Alternately, 5*3 = -D3 = 
Z/3Z X Z/2Z, so 5*3 could potentially arise from SSB in chiral p-wave superconductors. 

0(2) = f/(l) X Z/2Z broken down to S3 =< R120, flip > where denotes a counter- 
clockwise rotation by an angle u. 

In addition, hexagonal optical lattices can be manufactured using cold atoms. L. Duan, 
E. Demler, and M. Lukin proposed an implementation of Kitaev's code on the honeycomb 
lattice. [T7] Finally, a method for creating 5*3 anyons in Josephson junction arrays has been 
suggested. [IH] 
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2. As 



The conjugacy classes of are given by (), (12) (34), (123), (12345), (12354). The 
centrahzers of a representative element of each of these conjugacy classes are given by 
Z/2Z X Z/2Z (= group generated by (12)(34) and (13)(24)), Z/3Z (= group generated by 
(123)), Z/5Z (= group generated by (12345)), and Z/5Z (= group generated by (12354)) 
respectively. The dimensions of the internal Hilbert spaces of these fluxes arc given by the 
number of elements in the conjugacy class (a basis is given by linear combinations of elements 
in this conjugacy class). These dimensions are thus: 1, 15, 20, 12, 12 respectively. (The 
order of the conjugacy class times the order of the centralizer equals order of the group, as 
it must by elementary group theory.) 

First consider the pure charge sector. These are elements that arise from irreducible 
representations of A5. Since there are 5 conjugacy classes, there arc 5 irreducible represen- 
tations of A5. The dimensions of the irreducible representations are 1, 3, 3, 4, 5. Note that 
the sum of the squares of the dimensions of these representations is equal to l^lsl =60 as it 
should. 

The character table for ^5 is thus: 



Rep 





(12) (34) 


(123) 


(12345) 


(12354) 


Tl 


1 


1 


1 


1 


1 


T2 


3 


-1 





A 


*A 


T3 


3 


-1 





*A 


A 


74 


4 





1 


-1 


-1 


T5 


5 


1 


-1 









where A = —e^ — 

Now consider the dyons (elements of both flux and charge) - There are 4 possible fluxes 
other than the identity, labeled by the conjugacy classes of (12) (34), (123), (12345), and 
(12354) respectively. The possible electric charges for a particle with a flux corresponding 
to the conjugacy class of (12) (34) are given by representations of the centralizer - In this 
case, the centralizer is the Klein-4 Group: Z/2Z x Z/2Z, 
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Rep 


c 


(12)(34) 


(13)(24) 


(14) (23) 


Vi 

1 J- 


1 


1 


1 


1 




1 


-1 


-1 


1 




1 


-1 


1 


-1 


V4 


1 


1 


-1 


-1 



The possible electric charges associated with the flux corresponding to the conjugacy 
class of (123) are given by representations of the centralizer Z/3Z of this flux. Explicitly: 



Rep 


e 


(123) 


(132) 


Pi 


1 


1 


1 


P2 


1 


g27ri/3 


g47ri/3 


P3 


1 


g47ri/3 


g27ri/3 



Finally, consider the fluxes corresponding to the 5 cycles, (12345) and (12354). The 
electric charges in both cases are given by representations of Z/5Z, so the representation 
table will be listed only once: 



Rep 


e 


(12345) 


(13524) 


(14253) 


(15432) 


Pi 


1 


1 


1 


1 


1 


P2 


1 


g27rj/5 


g47rj/5 


g67ri/5 


g87ri/5 


P3 


1 


g47rj/5 


g37ri/5 


g27ri/5 


g67ri/5 


P4 


1 


g67ri/5 


g27ri/5 


g87ri/5 


g47ri/5 


P5 


1 


g87ri/5 


g67ri/5 


g47ri/5 


g27ri/5 



The particles are thus given by: 



30 



Particle 


Flux 


Charge 


Spin 


Dim Charge Repn 


Dim Internal Hilbert Space 


A 


V/ 


Ti 

1 





1 


1 


B 




V/ 


To 





3 


3 


C 





T-? 

' o 





3 


3 


D 





T4 





4 


4 


E 





Tt; 

o 





5 


5 


F 


(12)(34) 







1 


15 


G 


(12)(34) 


V-? 





1 


15 


H 


(12)(34) 




1 

2 


1 


15 


I 


(12)(34) 


1^ 


1 

2 


1 


15 


J 


(123) 







1 


20 


K 


(123) 


P'? 


1 

3 


1 


20 


L 


(123) 


r o 


2 
3 


1 


20 


M 


(12345) 


f 1- 





1 


12 


N 


(12345) 


P9 


1 

5 


1 


12 





(12345) 


P-^ 


2 

5 


1 


12 


P 


(12345) 


Pa 


3 
5 


1 


12 


Q 


(12345) 


P^ 


4 
5 


1 


12 


R 


(12354) 


P^ 





1 


12 


S 


(12354) 


P2 


1 

5 


1 


12 


T 


(12354) 


P3 


2 
5 


1 


12 


U 


(12354) 


P4 


3 
5 


1 


12 


V 


(12354) 


P5 


4 
5 


1 


12 



In summary then, there are 22 different types of particles. As a check, since the algebra 
is semisimplc, the sum of the squares of the dimensions of the internal Hilbert spaces for 
the particles should equal |^5p, the order of the quantum double. Now, the dimension of 
the internal Hilbert space is given by the product of the number of elements in the specific 
conjugacy class with the dimension of the irreducible representation. Thus the dimension is 
1^ + 32 + 3^ + 42 + 5^ + 52 + 4 * (15 * 1)2 + 3 * (20 * 1)2 + 2 * 5 * (12 * 1)^ = 3600. This does 
indeed equal \A^\^ = 60^ — 3600, as desired. 
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C 





E 


F 


G 


H 




A 


3 


c 





E 


F 


G 


H 


a 




i-B-E 


□ -E 


C-+4>4E 


E(C4[>t^ 


G-Htl 


F4H4I 


F4G-I 


i 


C 


D-E 


A-CtE 


EfOtE 


EfCfEH-: 


G-H+l 


F4i4l 


F4G-I 


D 


□ 


C+C-E 


B4D-E 


A-B-C+Mi 


B-tC-ID-2E 




F-G-H4I 


F4^4m- 


E 


E 


B«;-ri)-E 


B*C4D-E 


B+C-fDf!E 


i-B-iT-ZDtiE 


;f+g-hh 


=4!G4H-I 


F4C42HH 


F 


F 


G-Hfl 


Gfm- 


F-G+Kt 


i=-t<;-n-i 


i-Dti=-JF-:-li+_->.'4N4(>t-ii4C+q+StT-U-V 


B4C4CHE4rttStj4-:4L-M-N-0-P4a-R-54y4y4V 


a4C4D-E4G4l4j4K-L4MKj4.a4P4Q.4RtS-T44j4V 


_G 


S 


F+H-l 


F-Hfl 


F-GtH-n 


f^ZG^■^-\ 


E- 


-:i-E-H-l-J- -L-'.'-V-:i-:( :.-:.-!-T-U4V 


A4l>-ZE42GT-JM;tL4V4ft404'=4C.4Jl.+S4T-U-V 


B4C-D4E4F4:+J4*4L-M-rj-C-P-:i.-P-;-"tJ4V 


Jl 


H 


F4G-I 


F-Gfl 


F-G+t*t! 




B- 


-: -E-;-i-j-i 4 ( :-: -;-~-U4V 


R4C4D-E+.?+l4J4K+1.4M4N4t>4P-Cl4P4^T4UtV 


i4DflE42H4,4K4_-" ■ : ■ ■ -'-Tf-^-V 


1 


1 


F+G+H 


F+G-H 


F-Gtm* 
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FIG. 5: Fusion Rules for 

3. Fusion Rules: 

The fusion rules for can be found in the following pages. Once again every particle is 
its own antiparticle. Note that, unlike in the 5*3 case, there are multiple ways of obtaining 
some of the products - For instance, E and F can fuse to produce E in two distinct ways. 

IX. MAGMA PROGRAM AND DISCUSSION OF RESULTS 

The tables computed in this paper were produced using the MAGMA program included 
in Appendix 1. The MAGMA program computes the particles (and their properties - e.g. 
spin), braiding, and fusion rules for an arbitrary finite residual gauge group. The program 
can also explicitly compute the entries of the S matrix. As part of the fusionprinter( Group) 
procedure, the MAGMA program computes the change of basis S matrix and implements the 
Verlinde formula to determine the fusion coefficients based on the S matrix for an arbitrary 
group. 

In addition, the MAGMA program can compute the probabilities of obtaining any given 
fusion product used in many different topological quantum computation schemes. More 
detail on this connection is provided in the next section. 
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FIG. 6: Fusion Rules for 

The two examples included here were chosen mainly due to their possible applications to 
universal quantum computation (as described in the next section). In particular, since 5*3 is 
the smallest non-abelian group, a method for achieving universal quantum computation in 
that case has the best chance of being engineered in the lab or found in nature. However, 
certain other groups are interesting from either a physics theory or implementation point of 
view. In Appendix 2, there have been included some more examples of the fusion rules for 
the particles derived for other symmetry group quantum doubles that may be of interest. 
For example, from the physics theory perspective, in the two theories discussed above, every 
particle is its own antiparticle. However, this is not always the case and examples are 
included to illustrate this point. Furthermore, other dihedral groups, for example, could be 
interesting from the implementation perspective as they arise naturally from spontaneous 
symmetry breaking in p-wave superconductors. 

A. Majorana Particles 

From the fusion tables derived for S3 and A^, every particle is its own antiparticle, i.e. the 
particles are all Majorana. Most commonly, photons are their own antiparticles. However, 
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FIG. 7: Fusion Rules for A5 

Majorana particles of noninteger spin are currently an active area of research in physics, even 
beyond condensed matter applications. In particular, neutrinos are conjectured by some to 
be Majorana and SUSY would demand a fermionic photino as the fermionic superpartner of 
the photon (among many other Majorana particles). Finally, one suggestion is that WIMPS 
(a dark matter candidate, possibly the dominant component of dark matter mass in our 
universe) are also Majorana. 

X. QUANTUM COMPUTATION SCHEMES 

A. Universal Quantum Computation 

Universal computation: any unitary operation can be approximated to arbitrary accuracy 

B. Schemes for Topological Quantum Computation 

States used in a topological quantum computation scheme must be selected according to 
their fusion rules since the environment can distinguish between states fusing to different 
superselection sectors, causing decoherence. In addition, most implementation proposals 
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FIG. 8: Fusion Rules for A5 

use particle-antiparticle pairs as the fundamental building block since these interact trivially 
with other pairs when kept far apart from each other. Braiding computations are achieved 
by transporting a given pair in between other pair(s). 

If anyons from different superselection sectors are used in the actual quantum compu- 
tation scheme, then two distance scales come into play. In particular, the anyons must be 
kept far enough apart from each other that they cannot influence each other by local inter- 
actions (and thus superselection sectors are exhibited). Similarly, the environment must be 
far enough removed that it can only detect the cumulative long range quantum numbers of 
the whole system. (This provides protection against decoherence.) 

The prototype of anyon computation based on superselection sectors is the Fibonacci 
(also known as Yang-Lee) anyons, in which there are only two superselection sectors and 
1, with nontrivial fusion rule 1x1 = + 1. The dimension of the Hilbert space fusing to the 
vacuum increases exponentially in the golden ratio = — "^^^ associated R-matrix 
and its F-matrix conjugate FRF~^ then form a dense subgroup of SU{2). Thus any unitary 
operation can be approximated with arbitrary accuracy and universal quantum computation 
is possible. The Fibonacci anyons have spin given by | (This can be obtained by, for instance, 
solving the pentagon and hexagons explicitly for the R matrix and then using the generalized 
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d.AdBF(AB C) = 




FIG. 9: Fusion Probabilities. 

spin-statistics connection to infer the spins.) Thus this model is equivalent to SO{3)3 (the 
restriction of SU{2)3). 

1. Fusion Probabilities 

Finally, the probability results obtained for Fibonacci anyons in [19] can be extended to 
anyons based on the quantum groups of finite groups as well. Adopting Preskill's diagram 
convention of weighting by a compensating factor of \/(1a for each anyon-antianyon kink in 
a particle A's worldline, the dimension becomes equal to that of a topologically equivalent 
worldline. (Ordinarily every anyon-antianyon worldline zigzag would generate a factor of 
— where dA denotes the quantum dimension of particle A.) 

"A 

Consider the topological quantum computation schemes in which particle-antiparticle 
pairs are created from the vacuum. (This is standard in Kitaev's, Preskill's and Mochon's 
schemes.) Using the above diagram convention, Figures [9] and [10] show that the following 
formula for the probability that a particular fusion product occurs holds: 

FiAB ^C) = 

a Ads 
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dAdB¥(AB C] = ... = 




FIG. 10: Fusion Probabilities Cont'd. 
C. Possible Realizations of Topological Quantum Computation 

1. A5 

Ogburn and Preskill showed that = PSL{2,¥4) is sufficient for universal quantum 
computation under the assumption that a supply of magnetic ancillas is available (but not 
necessarily electric ancillas). [20] This sufficiency condition was subsequently generalized to 
any nonsolvable finite group. [21] Note that geometrically is the subgroup of rotations 
of the icosahedral group (the symmetries of a regular icosahedron - one of the five platonic 
solids), which in turn is a subgroup of 5*0(3). Out of the class of nonsolvable finite groups, 
A^ has the best chance of being engineered in the lab since it has the smallest order (60 
elements). It is unlikely to find a symmetry group with such high order in nature in a useful 
form. 

2. S3 

Since the smallest nonsolvable group is A^ and 1 = 60, it is desirable to analyze 
smaller groups. 6*3 = is the smallest nonabelian group {\Ss\ = 6 = ^), so in that 
sense a quantum computation scheme based on S3 would be optimal. S3 symmetry groups 
arise naturally in two dimensional lattices as the symmetry group of a hexagonal lattice 
(for example, graphene) so it is perhaps more reasonable to expect to find an 5*3 system 
conducive to topological gauge theories in nature. In addition, proposals have been suggested 
for engineering S3 in optical lattices or ion traps. 

Since there are 3-particle nonabelian subsystems of S3, one might hope that one would 
lead to universal quantum computation in a manner similar to the Fibonacci anyons. Un- 
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fortunately, the subsystems do not generate an SU(2)-dense 3 strand representation of the 
braid group as they do in the case of the Fibonacci anyons. In particular, consider a 3- 
particle subsystem of S3. Note that if we could also get rid of the B particle, then we obtain 
a subsystem with fusion rules given by the Fibonacci anyons. Unfortunately, the spins, and 
thus the R matrices, are not the same as in the Fibonacci case. In general, the SU{2)k 

1 / 27rz 

particles have spins that are specific multiples of (phase factors of e*=+2). 

For the moment, let P denote the particle with the nontrivial fusion rule (in all cases, 
the dimension of the internal Hilbert space of P is 2). N^p, where nP represents n P par- 
ticles, satisfies the recursion N^p = 2* A^^_2P + ^n-ip- The solution to this recursion is 
(1 + v^)" — (1 — V^)" 

N^p = . Thus the resulting subsystem has an exponentially increas- 

"'"^ 2^(2) 

ing Hilbert space dimension. In the case of P = G oi P = H, the particles have spin | or 
I and the fusion rules correspond to those of S0{3)i or its parity image. Unfortunately, 
50(3)4 is the restriction of SU (2)4 and it has been shown that SU (2)4 does not yield univer- 
sal quantum computation. [22j P = C or P = F have spin 0, so do not even have nontrivial 
R matrices. Thus the braid group representation described for the Fibonacci anyons would 
be trivial when considering combinations of F anyons. 

A method for universal quantum computation based on qutrits (as opposed to qubits) in 
5*3 that uses more than 3 particle types (including D fluxes and electric charges in addition 
to the vacuum A) was provided by Kitaev. The construction was generalized to solvable 
but non-nilpotent groups. [23] The major drawback of the method relative to the case is 
that it requires creation and fusion of electric charges, which may be more challenging than 
magnetic fluxes to create. In particular, the computational power of just conjugation alone 
will not give a Toffoli gate in the case of 5*3 as it does in the case of A^. 



3. The Dihedral Groups and p + ip Superconductors 

From the field theory point of view, time reversal or parity symmetry is broken in chiral 
p-wave superconductors so the electromagnetic U{1) group becomes minimally non-Abelian 
- U{1) X Z/2Z. Superconductivity results when this gauge group is spontaneously broken, 
and could give rise to non-Abelian anyons as, for instance, Dn = Z/nZ x Z/2Z. 
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XI. CONCLUSIONS 



In this paper, (2+l)-dimensional topological gauge theories based on finite gauge theories 
are analyzed with a view towards physical implementations and quantum computation. 
A program to compute the particles, their properties, the S matrix, and the fusion rules 
for a general finite group was written. The fusion rules for the quantum doubles most 
physically or computationally relevant were tabulated and the resulting tables were analyzed. 
In particular, the particles, fusion rules, and subsystems for two groups thought to be 
sufficient for universal quantum computation under certain circumstances - and S3 - 
were analyzed in some detail. The S3 theory with 8 particles was found to have many three 
particle subsystems. Unfortunately, these either led to trivial R matrices or a system given 
by 30(3)4, which has been shown to not lead to universal quantum computation since it 
does not produce a dense braid group representation. In addition, the particles, properties, 
and fusion rules for were tabulated and analyzed. 
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Appendix A: Appendix 

function ComputeSMatrix(A, B, alpha, beta, Gp) 
CC:=ConjugacyClasses(Gp) ; 

Cent := [Centralizer(Gp,Gp!CC[i] [3] ) : i in [l..#CC]]; 
total :=0; 

for j in [1 .. #Right Transversal (Gp, Cent [B] )] do 
for i in [1 .. #Right Transversal (Gp, Cent [A] )] do 
cj := RightTransversal(Gp, Cent [B] ) [j] ; 
ci : = Right Transversal (Gp , Cent [A] ) [i] ; 

cj :=cj"-l; 
ci:=ci~-l; 

valpha := ci"-l * cj * CC[B] [3] * cj"-l * ci; 
valbeta := cj''-l * ci * CC[A] [3] * ci^-l * c j ; 

if (Gp!cj * CC[B][3] * cj "-1)^-1 * (Gp! ci * CC[A] [3] * ci^-l)~-l 

*(Gp!cj * CC[B][3] * cj^-1) * (Gp! ci * CC[A] [3] * ci~-l) eq Id(Gp) 

then total :=total+ComplexConjugate((CharacterT able (Cent [A] ) [alpha] ) (valpha))* 
ComplexConjugate( (CharacterTable(Cent [B] ) [beta] ) (valbeta) ) ; 

end if ; 

end for; 

end for; 

total :=total/#(Gp) ; 
return(total) ; 
end function; 

function FusionCoef f icientsviaSMatrix(A, alpha, B, beta, C, gamma, Gp) 
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CC:=ConjugacyClasses(Gp) ; 

Cent := [Centralizer(Gp,Gp!CC[i] [3] ) : i in [l..#CC]]; 
fusionsum:=0; 

//sum over all centralizers D and representations of those centralizers delta 

for D in [1 . .#ConjugacyClasses(Gp)] , delta in [1 . .#ConjugacyClasses(Cent [D] )] do 

//since the number of conjugacy classes = number of representations 

//want delta up to # of repns of the centralizer corresponding to repn D Cent [D] 

if (ComputeSMatrix(l,D,l, delta, Gp) ne 0) 

then 

fusionsum:=fusionsum+ComputeSMatrix(A, D, alpha, delta, Gp) * 
ComputeSMatrix(B, D, beta, delta, Gp) * 

ComplexConjugate(ComputeSMatrix(C, D, gamma, delta, Gp))/ 

ComputeSMatrix(l, D, 1, delta, Gp) ; 
end if; 
end for; 

return fusionsum; 
end function; 

function FusionProb(A, alpha, B, beta, C, gamma, Gp) 
CC:=ConjugacyClasses(Gp) ; 

DA : =CharacterTable (Centralizer (Gp , Gp ! CC [A] [3] ) ) [alpha] [1] ; 
dimA:= (CC[A] [2])*DA; 

DB : =CharacterTable (Centralizer (Gp , Gp ! CC [B] [3] ) ) [beta] [1] ; 
dimB:= (CC[B] [2])*DB; 

DC : =CharacterTable (Centralizer (Gp , Gp ! CC [C] [3] ) ) [gamma] [1] ; 
dimC:= (CC[C] [2])*DC; 

fus:=FusionCoeff icientsviaSMatrix(A, alpha, B, beta, C, gamma, Gp) ; 
prob:= ((fus) * (dimC) ) / (dimA *dimB) ; 

return (prob) ; 
end function; 

function f usionprinter (Gp) 
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CC:=ConjugacyClasses(Gp) ; 

Cent := [Centralizer(Gp,Gp!CC[i] [3] ) : i in [l..#CC]]; 

for i in [1 . .#ConjugacyClasses(Gp)] , 

il in [1 . .#ConjugacyClasses(Cent [i] )] do 

printf " \n"; 

for j in [1. .#ConjugacyClasses(Gp)] , 

jl in [1 . .#ConjugacyClasses(Cent [j] )] do 

printf 

firstprint := 1; 
particlenum := 0; 

for k in [1. .#ConjugacyClasses(Gp)] , kl in 
[1 . .#ConjugacyClasses(Cent [k] )] do 
particlenum := particlenum + 1; 

if (FusionCoef f icientsviaSMatrix(i , il, j, jl, k, kl,Gp) ne 0) then 
if firstprint eq then printf "+"; 
end if ; 

firstprint := 0; 

if (FusionCoef f icientsviaSMatrixd , il, j, jl, k, kl,Gp) ge 2) then 
printf "\%o.", FusionCoef f icientsviaSMatrixd , il, j, jl, k, kl,Gp); 
end if; 

printf "\°/oO", particlenum; 
end if ; 
end for; 
end for; 
end for; 

return 0; 
end function; 

function spinprinter (Gp, A, alpha) //returns spin mod 1 
CC:=ConjugacyClasses(Gp) ; 
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CT : =CharacterTable (Centralizer (Gp , Gp ! CC [A] [3] ) ) [alpha] ; 

C := ComplexFieldO ; 

z: = (C!CT(CC[A] [3] ))/(C!CT[l] ) ; 

sl:=Log(z)/(2*Pi(C)*C.l) ; 

s : = (sl+ComplexConjugate (si) ) /2 ; 

return (s) ; 

end function; 

Appendix B: Other Interesting Quantum Doubles 

While 5*3 and Ar, were emphasized in the main body of the text, the same procedure was 
run on a number of other interesting non-abehan groups. The resulting particle tables and 
fusion rules for A4 are included in this appendix to demonstrate that (2+l)D topological 
gauge theories in which non-Majorana particles (particles that are not their own antiparti- 
cles) do indeed exist. S4 is included as the next highest order symmetric group. Finally, the 
second lowest order dihedral group D4 will be considered since dihedral groups could arise 
in some specific physical theories. 

1. A4 

A4 is the next smallest alternating group. Representatives for the conjugacy classes 
are given by: Id, (12) (34) (3 elements in the conjugacy class), (123) (4 elements in the 
conjugacy class), and (132) (4 elements in the conjugacy class) - Their centrahzers are 

A4 =< (12)(34),(234) >, Z/2Z x Z/2Z =< (13)(24), (12)(34) >, Z/3Z =< (123) >, and 
Z/3Z =< (123) >. The character tables of the centrahzers were all given in the text for 
another group, except of course A4 itself. The character table for A4 is: 
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Representations of A4 



Repn 


(1) 


(12)(34) 


(123) 


(132) 


Ki 


1 


1 


1 


1 


K2 


1 


1 


-K 


c 


«3 


1 


1 


C 


-1-C 


K4 


3 


-1 









where (2 + ^ + 1^0. 

Particles are specified by flux (a choice of conjugacy class) and charge (specifled by an 
irreducible representation of the centralizer of a representative element of the flux conjugacy 
class). Thus there are 14 particles, given by: 



Particle 


Flux 


Charge 


Spin 


Repn Dim 


Dim Internal Hilbert Space 


A 










1 


1 


B 










1 


1 


C 





f^s 





1 


1 


D 










3 


3 


E 


(12) (34) 


Vi 





1 


3 


F 


(12)(34) 







1 


3 


G 


(12)(34) 


V3 


1 
2 


1 


3 


H 


(12)(34) 


Va 


1 
2 


1 


3 


1 


(123) 


Pi 





1 


4 


J 


(123) 


P2 


1 
3 


1 


4 


K 


(123) 


P3 


2 
3 


1 


4 


L 


(132) 


Pi 





1 


4 


M 


(132) 


P2 


2 
3 


1 


4 


N 


(132) 


P3 


1 
3 


1 


4 



Note that the sum of the squares of the internal Hilbert spaces is 3*1^+32+4*32+6*42 = 
144 = 1^4]^, the dimension of the quantum double. 

Note that B and C are antiparticles, so this is the flrst example we have of a theory in 
which not every particle is its own antiparticle. 
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A 


1 




J5 


I 


i 


s 


Ji 


i 


i 




1 


M 


Ji 


A 


A 


3 


C 


D 


E 


F 


G 


H 


1 


J 


K 


L 


M 


N 




B 


C 


A 





E 


F 


G 


H 


K 


1 


J 


N 


L 


M 


X 


r 

y_ 


A 


B 


D 


E 


F 


G 


H 


J 


K 


1 


M 


N 


L 


D 


D 


D 


D 


A+B+C+2D 


F+G-H 


E+G+H 


E+F+H 


E+F+G 


l+J-K 


l+J+K 


l+J+K 


L+M+N 


L-M+N 


L+M+N 


J 


E 


E 


E 


F-fGfH 


A+BfC+2E 


CHG+H 


[HF+H 


[HF+G 


l+J+K 


l+J+K 


i+J-K 


L+M+N 


L+M+N 


L+M+N 


i 


F 


F 


F 


E+C^H 


D+G+H 


A+E+C+2F 


D+E+H 


D+E+G 


l+J+K 


l+j+K 


l+J-K 


L+M+N 


L+M+N 


L+M+N 


G 


G 


G 


G 


E+F+H 


D+F+H 


CHE+H 


A+fl+<:+2G 


[>+E+f 


l+J+K 


l+J+K 


l+J+K 


L+M+tJ 


L+M+N 


L+M+N 


H 


H 


H 


H 


EfF+G 


D+-F+G 


D+E+G 


D+E+F 


A+B+C+2H 


l+J-K 


l+J-K 


l+J+K 


L+M+« 


L+M+N 


L+M+N 


\ 


1 


K 


J 


IfJ+K 


liJ+IC 


l+J+K 


l+J+K 


l+J+K 


2L+M+N 


L+2M-N 


L+M+2N 


A+D+E+F+G+H 


C+&+E+F+G-H 


B+D-E+F+G+H 


J 


J 


1 


K 


IfJ+K 


l+JiK 


l+J+K 


l+J+K 


l+J+K 


L+2M+N 


L+M+2N 


2L+M+N 


C+D+E+F+G+H 


a+[>+E+F+^-H 


A-D-E+F+G+H 


K 


K 


J 


1 


l+J+K 


l*J+K 


l+J+K 


l+J+K 


l+J+K 


L+M+2N 


2L+MtN 


L+2M+N 


&+D+€+F+&+H 


A+CHE+F+G+H 


C+D-E+F+G+H 


I 


L 


N 


M 


L+M+N 


L+M+14 


L+M+N 


L+M+N 


L+M+N 


A+[}+E+F+G+H 


C+D+E+F+G+H 


B+D+E+F+G+H 


21+J+K 


I+2J+K 


I+J+2K 


M 




L 


N 


L+M+N 


L+M+N 


L+M+N 


L+M+N 


L+M+W 


C+D+E+F+Gr+H 


B+D+E+F+G+H 


A+[>+E+F+<5+H 


l+2JtK 


I+J+2K 


2ltJ+K 


H 


H 


M 


L 


L+M+N 


L+M+N 


L+M+N 


L+M+N 


L+M+W 


B+D+E+F+GH-H 


A+D+E+F+G+H 


C+D+E+F+&+H 


I+J+2K 


21+J+K 


I+2J+K 



FIG. 11: Fusion Rules for 

2. 54 

The next highest order symmetric group is 5*4, with 4! = 120 elements. There are five 
conjugacy classes (specified by the cycle decomposition) - Id, (12), (12) (34), (123), and 
(1234). The centralizers for these conjugacy classes are, respectively, 5*4, D4, TLjTL x TLjlL^ 
Z/3Z, and Tj/ATj. The character tables for all of these groups (with the exception of 6*4 and 
D4) have been given in the text, so will not be repeated here (the same representation labels 
will be used in the enumeration of the particles). 

The character table for 6*4 is: 



46 



Irrep 




V/ 


(13)(24) 


(12) 


(123) 


(1234) 


A. -L 


1 


1 


1 


1 


1 


Y2 


1 


1 


-1 


1 


-1 


X3 


2 


2 





-1 





X4 


3 


-1 


-1 





1 


X5 


3 


-1 


1 





-1 



The character table of D4 is: 



Irrep 





(13)(24) 


(14)(23) 


(24) 


(1234) 


Ti 


1 


1 


1 


1 


1 


T2 


1 


1 


1 


-1 


-1 


T3 


1 


1 


-1 


-1 


1 


T4 


1 


1 


-1 


1 


-1 


T5 


2 


-2 












There are 21 particles. 
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Particle 


Flux 


Charge 


Spin 


A 





A, 1 





B 


V/ 







C 










D 


V/ 







E 


V/ 


AO 





F 


(12)(34) 


Ti 





G 


(12)(34) 







H 


(12) (34) 








I 


(12) (34) 


T4 





J 


(12) (34) 


T5 


1 

2 


K 


(12) 


/-L 





L 


(12) 







M 


(12) 


to 


1 

2 


N 


(12) 


Va 


1 

2 





(123) 


1 J- 





P 


(123) 


Pi 


1 

3 


Q 


(123) 


P-i 
ro 


2 
3 


R 


(1234) 


Pi 





S 


(1234) 


P2 


1 

4 


T 


(1234) 


h 


1 

2 


U 


(1234) 


pi 


3 
4 



From the fusion rule table, every particle is its own antiparticle. 
3. L>4 

The conjugacy classes are given by Id, (13)(24) (1 element in the conjugacy class), (14)(23) 
(2 elements in the conjugacy class), (24) (2 elements in the conjugacy class), (1234) (2 
elements in the conjugacy class) - The centralizers are respectively D4, D4, Z/2Z x Z/2Z =< 
(12)(34), (24) >, Z/2Z x Z/2Z =< (13), (24) >, Z/4Z =< (1234) >. 

Thus there are 22 particles. 
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FIG. 12: Fusion Rules for 
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FIG. 13: Fusion Rules for ^4 
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FIG. 14: Fusion Rules for 
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FIG. 15: Fusion Rules for ^4 
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As a check, note that the sum of the dimensions of the internal Hilbert spaces of the 
particles is 8 * 1^ + 14 * 2^ = 64 = the dimension of the quantum double. 

Note by inspection that each particle is again its own antiparticle in this theory. 

a. Possible Physical Realization 

Note that D4 = Z/4Z k Z/2Z, so one possible physical reahzation could arise from SSB 
in chiral p-wave superconductors. 
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FIG. 16: Fusion Rules for D4 
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FIG. 17: Fusion Rules for D4 
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